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Abstract 



This paper introduces the FLOSS Free Libre Open Source 
Software [VEMSA3D] , a contraction of "Visual Electromagnetic 
Simulator for 3D Antennas", which are geometrically modeled, 
either exactly or approximately, as thin wire polygonal stru- 
ctures; presents its GUI Graphical User Interface capabili- 
ties, in interactive mode and/or in handling suitable formed 
antenna data files; demonstrates the effectiveness of its use 
in a number of practical antenna applications, with direct 
comparison to experimental measurements and other freeware 
results; and provides the inexperienced user with a specific 
list of instructions to successfully build the given source 
code by using only freely available IDE Integrated Develop- 
ment Environment tools-including a cross-platform one. The 
unrestricted access to source code, beyond the ability for 
immediate software improvement, offers to independent users 
and volunteer groups an expandable, in any way, visual an- 
tenna simulator, for a genuine research and development work 
in the field of antennas, adaptable to their needs. 

Keywords 



FLOSS, antennas, modeling, simulation 



Introduction 

A lot of amazing visual EM 
software simulators, both com- 
mercial and freeware, exist 
for many years now. However, 
to the best of authors' know- 
ledge, none of these simula- 
tors has been ever released 
under an Open Source license- 
paid, free gratis or free li- 
bre. This situation, as it 
was discussed recently, pre- 



cludes the independent users 
and volunteer groups with 
limited resources, like the 
authors' nonprofit group, from 
the genuine, state-of-the- 
art, scientific research, sin- 
ce, neither the code improve- 
ment, nor its adaptation to 
specific needs, is possible 
[ 1 ], [ 2 ]. 

For that reason, the au- 
thors decided to develop and 
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release under the GNU public 
license their own Visual Elec- 
tromagnetic Simulator for 3D 
Antennas [VEMSA3D] , although 
it certainly is a less well- 
equipped software application 
for the moment-but still a 
fully expandable one. That re- 
lease is now perfectly per- 
missible, because in 2005 
NASA released in the public 
domain the FORTRAN source co- 
de of the well-known MoM Me- 
thod of Moments Thin-Wire Com- 
puter Program, by J. H. Ri- 
chmond [3], [4]. This is ex- 
actly the code on which the 
authors' group based its free- 
ware simulators: [RichWire], 
a CLI Command Line Interface, 
and [DA], a MS Quick-Win FOR- 
TRAN derivative of it. These 
two simulators are under un- 
interrupted revision, impro- 
vement, expansion, and rede- 
velopment, since a long time 
ago [1], [2]. Therefore, the 
[Rich-Wire] FORTRAN code was 
translated-entirely, line-by- 
line, without using any paid 
or free translator-to C++, to 
form the core of scientific 
EM computations in [VEMSA3D] . 

On the other hand, the au- 
thors' group requirements for 
scientific accuracy in visual 
representation of the produ- 
ced EM simulation antenna re- 
sults from its simulators, we- 
re already enforced the soft- 
ware development of the Vir- 
tual Antennas, that is the 
Virtual Antennas laboratories. 



in VRML [5], [6], the FLOSS 

application [RadPat4W] for an- 
tenna radiation pattern pre- 
sentation, in MS VB6 [2], 
[7], as well as, the recently 
developed visual antenna ap- 
plication for the Wolfram De- 
monstrations Project in Ma- 
thematica [8] . The visualiza- 
tion ideas implemented in the 
aforementioned graphics ap- 
plications were also expres- 
sed-from the scratch-in C++, 
to use the Open Source cross- 
platform [wxWidgets] library 
with OpenGL and form the core 
of scientific EM graphics in 
[VEMSA3D] [9] . 

The authors, having taken 
into account that no famil- 
iarization with software use 
is possible without getting 
hands-on experience, restrict 
themselves to a brief discus- 
sion of the current [VEMSA3D] 
characteristics. The code, the 
antenna applications data, as 
well as, any other informa- 
tion, referenced open bugs or 
future code releases, will be 
always available in authors' 
group repositories, at 
"www. antennas . gr /floss" 
and in GoogleCode website, at 
"http : //code . google . comp/p/rg 
a/" . 

GUI Interactive Mode of 

Operation 

It is assumed that the 
user has already some experi- 
ence in the sketching of a 
polygonal wire outline model, 
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for an antenna under consid- 
eration, consisting of num- 
bered wire segments and no- 
des, with their 3D space co- 
ordinates, including the po- 
sitions of any antenna cir- 
cuit elements, that is volt- 
age generators and lumped 
loads, and s/he wants then to 
key-in these model and cir- 
cuit data into [VEMSA3D] us- 
ing its GUI Graphical User 
Interface in interactive mo- 
de . 

The GUI main window is 
shown in Fig. 1, while Fig. 2 
shows the menu items along 
with their available submenu 
options numbered from 1 to 7. 

The GUI itself is divided 
in three panels named [An- 
tenna in space (3D)], [Anenna 
Elements], and [Data-sets]. 



The function of each panel is 
briefly described in the fol- 
lowing . 

The [Antenna in space 
(3D)] viewing panel is used 
to project all the generated 
3D and 2D graphics. By de- 
fault, the application starts 
with a simple linear dipole 
antenna in space, as it is 
shown in Fig. 1 from a view- 
point on the diagonal of tri- 
hedral angle OXYZ, where the 
usual letters associated with 
the Cartesian axes do not ex- 
ist. Instead, a one-to-one 
correspondence implicitly ex- 
ists between (X, Y, Z) axes 
and (R, G, B) colors [5]. The 
3D image of any antenna can 
be manipulated through the 
[Study] options, as shown in 
Fig. 2.3. 




Fig. 1: GUI: The main window 



TUESDAY 23 DECEMBER 2014 



Vl-9 



FUNKTECHNIKPLUS # JOURNAL 



C.A. KOUTSOS, N.I. YANNOPOULOU, P.E. ZIMOURTOPOULOS 



| View Study Graph Simulation Options Help 



New 

Open... 

Save 

Save As... 

Import Geometry... 
Export Geometry... 



Study Graph Simulation 
Tool Bar 2 
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Help 

Viewing Options... 0 
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Options MIS 
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Simulation Options Help 
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Antenna segments in space (3D) 

E Radiation Pattern in space (3D) 
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Fig. 2: Main window: The unfolded menu items 



The [Antenna Elements] pa- 
nel contains four tabs, cor- 
responding to model and cir- 
cuit data, in four self-ex- 
plained building blocks. In- 
stead of attempting to de- 
scribe the block data in gen- 
eral the presentation conti- 
nues with the default data: 

[Nodes] : 2, with coordi- 
nates 1: (0.2, 0.1, 0.25), 2: 
(0.2, 0.1, -0.25), in meters 
[m] or in terms of wavelength 
A [wl]. 

[Wires] : 1, starting at 
node: 1, ending at node: 2, 
that is a directed segment. 

[Segments] is either: (a) 
a positive integer, (b) a 
zero or (c) a number between 
0.00 and 0.25, to respec- 
tively define the division of 
this wire in segments: (a) in 
the indicated number of seg- 
ments, (b) in a calculated 
number of segments of length 
no longer than 0.05A (that is 
10 segments, in this case), 



(c) in segments with length 
no more to that number. No- 
tably, no segment can be 
longer than 0.25A [3]. 

[Generators] : 1, connected 
in the wire:l, at a distance 
of 0.5 times this wire 
length, from its starting no- 
de, with rms value of 1^0° 
[V] . There can be only one 
generator on each wire. 

[Loads]: Has the same stru- 
cture as [Generators] . Resis- 
tors, inductors and capaci- 
tors are inputted in [Q] , [H] 
and [F], respectively. It is 
empty in this case, since no 
lamped loads exist. 

Notably, the GUI has been 
developed in such a way that 
efficiently supports the in- 
teractive handling of the an- 
tenna modeling by the user on 
the screen. Thus, the selec- 
tion of an element in any of 
these four tabs simultane- 
ously highlights this element 
on the viewing [Antenna in 
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space (3D)] panel-if the [Se- 
lect] option under [Study] 
(Fig. 2.3) has been already 
chosen by the user, and this 
is the default selection in 
any case-and conversely: the 
selection of any element on 
this viewing [Antenna in spa- 
ce (3D)] panel highlights the 
element in its table. Conse- 
quently, a full control on 
the antenna modeling struc- 
ture is achieved. 

The [Datasets] panel may 
contain independent sets of 
model and circuit data for va- 
rious antennas. Such a Data- 
set has its own [Dataset Si- 
mulation Options], which are 
chosen through submenu of 
Fig. 2.6, as the default se- 
lection of them is shown on 
the left part of Fig. 3. 

A [Wire conductivity] of -1 



corresponds to perfect wire 
conductivity, while the [Spa- 
ce conductivity] of 0 and the 
[Space dielectric constant] 
of 1 defines the Free Space 
EM environment. The [Integra- 
tion steps] accepts a posi- 
tive number for the steps of 
approximated integration used 
in the MoM impedance computa- 
tions-a zero means closed- 
form integration. 

Notably, the [Maximum seg- 
ment length] affects the di- 
vision of wires: the smaller 
this number is, the more seg- 
ments will be considered. It 
is currently known that there 
is an implied, still program- 
matically unimproved, priori- 
ty of the number of segments 
defined in [Wires] tab, over 
this selection-that is a bug, 
which may crash the applica- 
tion . 



Dataset Simulation Opti 


ons X 


Operating frequency : 


300 


[MHz] 


Wire radius : 


0.001 


[m] 


Wire conductivity : 


1 _1 


[MS/m] 


Speed of light : 


3e4008 


[m/sec] 


Space conductivity : 


0 


[S/m] 


Space dielectric constant : 


1 


]B 


Maximum segment length : 


0.05 


Pambda] 


Integration steps : 


0 


]0 


3D angle step : 


5 


[deg.] 


2D angle step : 


1 


[deg.] 



OK ~j | Cancel | | Default ] 



Status 



20:20:35.453 


Simulation started. 


20:20:35.453 


Running Dataset : Dataset 


20:20:35.687 


SEG started. 


20:20:35.703 


SEG finished. 


20:20:35.703 


SORT started. 


20:20:35.750 


SORT finished. 


20:20:35.750 


SGANT started. 


20:25:32.718 


SGANT finished. 


20:25:32.718 


GANT1 started. 


21:03:11.906 


GANT1 finished. 


21:03:11.906 


PATTERNS started. 


21:03:22.500 


PATTERNS finished. 


21:03:22.500 


Simulation finished. 


21:03:22.500 


Total time: 00:42:47.047 



Close 



Fig. 3: Windows: [Dataset Simulation Options] - [Status] 
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The EM simulation of an an- 
tenna starts by selecting [Si- 
mulate Current Dataset], un- 
der [Simulation] (Fig. 2.5), 
while its progress is shown 
in the [Status] window, such 
as the one at the right of 
Fig. 3, which concerns the 
default dipole. Notably, mul- 
tiple Datasets can be simu- 
lated, one after the other, 
through the [Simulate Marked 
Datasets] shown in Fig. 2.5. 

The course of the simula- 
tion may be roughly described 
in five steps: 

(1) Wires are automatically 
segmented, using the provided 
parameters, to produce the fi- 
nal model structure of points 
and segments 

(2) Model structure is ana- 
lyzed and adjacent segments 
are combined to form dipole 
current modes with sinusoidal 



distribution 

(3) Impedances are calculated 
and the MoM algebraic system 
of equations is formed 

(4) The system of equations 
is solved and the segment cur- 
rents are calculated 

(5) 3D and 2D radiation pat- 
terns are calculated, as well 
as other useful results like 
those shown in Fig. 4. 

In Fig. 5, the [Viewing 
Options] under the [Options] 
(Fig. 2.6) is given with all 
of the available options for 
the graphics. Through [Graph] 
menu item (Fig. 2.4), a vari- 
ety of 3D and 2D plots, which 
include either the normalized 
electric far field E radia- 
tion pattern or the relative 
radiation intensity U in dB, 
as well as their 0 and cp 
parts, can be illustrated. 




Fig. 4: Window: [Results] - Table: [Power] 
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The [Antenna segments in 
space (3D)] illustrates the 
final antenna structure with 
points, segments and segment 
currents in amplitude and 
phase, as shown in Fig. 6, 
where the resulting current 
distribution on default di- 
pole it seems to be sinu- 
soidal, that is as it was ex- 
pected to be. The well-known 



2D intersections of the 3D 
radiation patterns by the 
three main planes yOz, xOy, 
zOx, as well as these 3D pat- 
terns themselves are also 
shown for the default dipole. 

Notably, modeled antennas 
can be imported and exported 
in [RichWire] data file for- 
mat (Fig 2.1) . 




Fig. 5: [Viewing Options] 




Fig. 6: GUI: Default dipole results 

Build under MS Windows XP 

To build [VESMA3D] MS Win- IDE CodeBlocks with GCC/Mingw 
dows executables, the use of is suggested [10] : 
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Set below as X:\ the system disk, e.g. X:\ -> C:\ 

1 Download (~20MB) , Install (~110MB) : 

http : //sourceforge . net/pro jects/codeblocks/f iles/Binaries/ 
[View all files] [Binaries] [8.02] 
codeblocks-8 . 02mingw-setup . exe 
Select: [Full], [Run], Compiler: [GCC/Mingw] 

2 Download wxWidgets (~20MB) , Extract (~120MB) : 

http : //sourceforge . net/pro jects/wxwindows/f iles/2 . 8 . 11/ 
[View all files] wxWidgets-2 . 8 . 11 . zip 
wxWidgets-2 . 8 . 11 . zip-> X:\ [Yes to All] 

3 IDE CodeBlocks, Open: 

X : \wxWidgets-2 . 8 . ll\include\wx\msw\setupO . h 
Replace, at line 1006: 

#def ine wxUSE_GLCANVAS 0 -> 1 
Repeat it, to file: setup. h 

4 At system variable [path] : 

[My Computer] [Properties] [Advanced] 

[Environment Variables] [System Variables] 

Path [Edit ] [Variable value], Add: 

" ; X : \Program Files\CodeBlocks\MinGW\bin; " 

5 Command Prompt : >cd\ 

>cd X : \wxWidgets-2 . 8 . ll\build\msw 
>mingw32-make -f makefile. gcc MONOLITHIC=0 
S HARED =0 UNICODES USE_OPENGL=l 
BUILD=release 

Re-Command, with option: 

BUILD=release -> BUILD=debug 

6 Download (~100KB) , Extract (~700KB) : 
http : //code . google . com/p/rga/ downloads/ 

VEMSA3D_source_l . zip -> X : \ 

7 IDE CodeBlocks, [File] [Open] : 

X : \VEMSA3D\build\win-cb\VEMSA3D . workspace 

* If system disk is not C:\ then: 

(#) [Project] [Build options. . .] [VEMSA3D] 

* [Search directories] Correct in each of: 

* [Compiler] , [Linker] , [Resource compiler] 

* C:\ -> X : \ [Yes] 

* Repeat from (#) , for [Release] 

* Repeat from (#) , for [Debug] 

[Build] [Select target] [Release] [Build] 

[Build] [Select target] [Debug] [Build] 
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This process results in 
[VEMSA3D] executables of ~35 
MB Debug and of ~5MB Release 
versions, which run under NT4, 
W2K and WXP, at least. 

Alternatively, the use of 
the freeware IDE MS Visual 
Studio Express with C++ 2008 
Compiler is also suggested. 
For that, first download MS 
Visual Studio 2008 Express 
iso-image from Microsoft web- 
site, burn it into a CD-R, 
and install it. Then, down- 
load from our repositories 
the file [VEMSA3Dfiles-4win- 
cb.zip], extract it, and fol- 
low the included setup in- 
structions. This process re- 
sults in [VEMSA3D] executa- 
bles of -5MB Debug and of 
~2MB Release versions, which 
run under W2K and WXP, at 
least, but definitely do not 
run under NT4 . 

Practical Antenna Applica- 
tions and Results 

This section presents the 
implementation of [VEMSA3D] to 
produce eight antenna models, 
from simple to more compli- 
cated structures. The input 
data were imported through 
[Import Geometry] of the [Fi- 
le] menu (Fig. 2.1). 

Fig. 7 shows an array of 2 
dipoles for operation at the 
frequency of 1111 MHz, dis- 
tanced by 0.85A, constructed 
by bare copper wire of 1 mm 
(0.0037A) radius, and mea- 
sured. In the same figure, 



the 3D radiation intensity 
pattern and its 2D main plane 
cuts are shown [11] . The con- 
tinuous line represents [VEM- 
SA3D] results, the dashed li- 
ne, analytically produced pat- 
terns, and the dots, measure- 
ments made using the authors' 
group VNA system [12], [13]. 

Fig. 8 illustrates the re- 
sults for a constructed and 
measured improved Hentenna mo- 
del, at 1110 MHz, with height 
of A/2, width of A/6, and 
with the active element at a 
distance 7A/12 from its bot- 
tom [14] . 

The third antenna consists 
of a A/4 monopole (5.83 cm) 
at 1286 MHz over a circular 
counterpoise of 14 cm radius. 
The monopole was constructed 
by 2 mm diameter copper wire 
and the counterpoise has been 
printed on a circuit board of 
29 cm x 29 cm, as a circle 
with four radials of 2 mm 
width [15]. In Fig. 9, the 
left 3D and 2D patterns are 
for vertical polarization, 
while the right ones for ho- 
rizontal. For a better repre- 
sentation, the patterns have 
been normalized. 

A somehow more complicated 
antenna is shown in Fig. 10. 
The left 3D and 2D patterns 
correspond to the vertical po- 
larization of a monopole over 
circular counterpoise with 16 
radials, while the right co- 
lumn patterns are the corre- 
sponding 3D and 2D patterns 
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for the same antenna, but with 
an additional grounded disc 
below it, and without elec- 
trical connection to the coun- 
terpoise antenna, at a dis- 
tance of 0.9 cm. The grounded 



disc was constructed from thin 
copper sheet and was modeled 
in [VEMSA3D] with 5 concen- 
tric circles and 64 radials 
[15]. 




Fig. 7: Array of 2 dipoles: 73 dipole current modes 




— VEMSA3D •••Measurements 



Fig. 8: Hentenna: 109 dipole current modes 
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Fig. 9: Monopole over a counterpoise with 4 radials: 
167 dipole current modes 




— VEMSA3D •••Measurements 



Fig. 10: Monopole over a counterpoise with 16 radials: 
419 dipole current modes. 

Monopole over a counterpoise with 16 radials over 
a grounded disc: 2594 dipole current modes 
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In Fig. 11, the results of 
[VEMSA3D] for a commercial 
UHF antenna are shown [16] . 
The antenna model is presen- 
ted in detail and separately, 
for the active element sys- 
tem, as well for the wire- 
frame reflector. The patterns 
correspond to the center fre- 
quency of its operation at 
650 MHz. 

Fig. 12 illustrates the 
model of a small jet airplane 
from the well-known freeware 
[4NEC2X] antenna simulator 
[17] . Simulation was carried 
out at the frequency of 10 
MHz with the same number of 
points and segments in both 
simulators. There is a good 
agreement between them for 



the radiation intensity pat- 
terns although some deviation 
in input impedance and direc- 
tivity is observed. 

Fig. 13 shows the results 
for a horn antenna with di- 
mensions proposed by K. Pitra 
and Z. Raida [18]. A small 
bow-tie feeder with a trian- 
gular perimeter of 0.47A and 
flare angle 37.50°, is con- 
sidered. The antenna was ini- 
tially simulated at the fre- 
quency of 40 GHz in both 
[4NEC2X] and [VEMSA3D] simu- 
lators, with 1205 segments to 
be consistent with the re- 
stricted maximum number of 
1500 segments of [4NEC2X] . 
There is a good agreement be- 
tween the produced results. 




Fig. 11: Typical commercial 
723 dipole current 



TV UHF antenna: 
modes . 
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[VEMSA3D] 



Fig. 12: Jet airplane modeled in [4NEC2X] and [VEMSA3D] : 
391 dipole current modes 



Finally, the most compli- 
cated antenna model, which is 
presented, corresponds to the 
same horn antenna at the same 
frequency of 40 GHz, with a 
dense wire-frame consisting 
of 3266 points and 4362 seg- 
ments. The model is shown in 
Fig. 14. The total number of 
dipole current modes is 5458, 



as shown in Fig. 15. The pro- 
cess time ranges from ~40 min 
in an AMD Phenom X2 550 3.11 
GHz CPU to -220 min, in an 
Intel Pentium 4 1.7 GHz CPU. 
Fig. 16 illustrates the re- 
sulting radiation patterns 
and Fig 17 the current ampli- 
tude and phase. 
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Ji 

w 

[4NEC2X] 




[VEMSA3D] 




Impedance | 


30 -i 


0.06 


S.W.R. 50 j 


1.67 




Frequency 


1 40000 


Mhz 


Wavelength 


| 7.e-3 


mtr 


Norm-All : 


: 6.53 dBi 






Resistance [Ohm] 

4 . 673E+001 

Reactance [Ohm] 

-3 . 929E+000 

VSWR [] 

1.1115 

Directivity 
6. 702E+000 [dBi] 



Fig. 13: Horn antenna with bow-tie feeder modeled in 
[4NEC2X] and [VEMSA3D] : 1370 dipole current modes 




Fig. 14: Horn antenna with bow-tie feeder: 4362 segments 
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Points Segments Dipole Current Modes ; input Sources Power 





Point 1 


Point 2 


Point 3 


Segment A 


Segment B 


Amplitude [A] 


Phase [deg.] 


5455 


1121 


3263 


1122 


4319 


4320 


1 . 053E-002 


-170.85 


5456 


1121 


3264 


1123 


4321 


4322 


1 . 042E-002 


-172.78 




1122 


3265 


1123 


4323 


4324 


1 . 615E-004 


-127.73 


^5458 


) 1121 


3266 


1118 


4325 


4326 


2 . 602E-002 


12.29 



Points Segments Dipole Current Modes 



Total Input Power (Real) 


2 . 543E-002 [W] 


Total Input Power (Imag) 


-5 . 540E-003 [W] 


Total Radiated Power 


2 . 543E-002 [W] 


Total Dissipated Power 


0 . 000E+000 [W] 


Power Efficiency 


100% 


Power Gain 


5 . 719E+000 [] 




7 . 573E+000 [dBi] 




5 . 424E+000 [dBd] 


Directivity 


5.719E+000 [] 




7 . 573E+000 [dBi] 




5 . 424E+000 [dBd] 



Fig. 15: Horn antenna with bow-tie feeder: 
5458 dipole current modes 




Fig. 16: Horn antenna with bow-tie feeder: 
Radiation intensity 
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Fig. 17: Horn antenna with bow-tie feeder: 
Current Amplitude and Phase 



Conclusion 

The first stable version 
of a visual EM simulator has 
been developed and released 
as FLOSS. The computational 
results of its use were found 
to be in good agreement with 
experimental measurements as 
well with the comparable free- 



ware simulator [4NEC2X] . Sin- 
ce the number of possible mo- 
difications and additions to 
the attributes of [VEMSA3D] 
seems to be endless, no at- 
tempt will be made to suggest 
a particular direction for its 
future development: any con- 
tribution from the antenna 
community is very welcomed. 
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A faded synthesis of an anthemion rooted in a meandros 

The thirteen-leaf is a symbol for a life tree leaf. 

"Herakles and Kerberos", ca. 530-500 BC, 
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Museum of Fine Arts, Boston. 
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The simple meandros is a symbol for eternal immortality. 

"Warrior with a phiale", ca. 480-460 BC, 

by Berliner Maler, 

Museo Archeologico Regionale "Antonio Salinas" di Palermo, 
commons .wikimedia . org/wiki/File : War rior_MAR_Palermo_N 12134. j pg 



